Phonon Hall effect in ionic crystals in the presence of static magnetic field 
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We study phonon Hall effect (PHE) for ionic crystals in the presence of static magnetic field. Using 
Green-Kubo formula, we present an exact calculation of thermal conductivity tensor by considering 
both positive and negative frequency phonons. Numerical results are shown for some lattices such 
as hexagonal lattices, triangular lattices, and square lattices. We find that the PHE occurs on 
the nonmagnetic ionic crystal NaCl, although the magnitude is very small which is due to the tiny 
charge-to-mass ratio of the ions. The off-diagonal thermal conductivity is finite for nonzero magnetic 
field and changes sign for high value of magnetic field at high temperature. We also found that the 
off-diagonal thermal conductivity diverges as ztf/T at low temperature. 
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I. INTRODUCTION 



The magnetic field dependence of heat conductiv- 
ity in paramagnetic dielectric crystals has been found 
experimentalljii. This effect is known as "phonon Hall 
effect" . Applying magnetic field in the direction normal 
to the temperature gradient, the authors have discov- 
ered the appearance of a transverse thermal current in 
Tb3Ga50i2 and the effect was confirmed in Ref. 0. Due 
to the Lorentz force on the electrons, electronic Hall ef- 
fect can be easily understood. However for phonons there 
is no such direct coupling with the magnetic field. The 
underlying mechanism that determines this effect is the 
spin-phonon interaction (SPI)^""— of phonons and param- 
agnetic ions Tb+^. This effect is an analog of the anoma- 
lous Hall effect (AHE) in paramagnetic phase. The the- 
ory of this phenomenon was proposed in Refs.0-H where 
the authors consider the spin-phonon interaction and 
solve the problem perturbatively. The same effect was 
also studied in four-terminal junctions using nonequilib- 
rium Green's function (NEGF) approach^. In Ref. [l^ 
the authors explained that molecules with rotary degree 
of freedom can also have the same effect as PHE through 
renormalization of acoustic waves. This effect for gases 
is knows as Senftleben-Beenakker effect. Very recently, 
the topological nature of PHE and an associated phase 
transition was found in the paramagnetic dielectrics^i. 

All the previous work on PUE^^^dL are based on the 
paramagnetic dielectrics; and the recent workJ^^ about 
the magnon Hall effect in the ferromagnetic insulator 
could also include the PHE to the thermal Hall effect. 
We know that all the materials in the previous experi- 
mental or theoretical work on PHE are paramagnetic or 
ferromagnetic ionic crystals; therefore, a question arises 
promptly: whether the PHE can occur in a general non- 
magnetic ionic crystal? From a theoretical point of view, 
it is highly desirable to study both the existence and 
properties of the PHE in general ionic crystals with an 
applied magnetic field, which could also guide us to do 
the further experimental study on PHE. 



As there is no magnetization effect in the nonmagnetic 
ionic crystal, intuitively one can not imagine the exis- 
tence of PHE in the sample. The Hamiltonian in the 
harmonic approximation for the ionic crystal in the pres- 
ence of a magnetic field is quadratic and is very similar 
to the form in Refs. Because of the quadratic nature 
of the problem it is possible to diagonalize the Hamilto- 
nian in terms of creation and annihilation operators. In 
Ref. d the authors consider the Hamiltonian which is not 
positive definite and omit the contribution of negative 
frequency to the phonon Hall conductivity which was in- 
correct. Here we carefully reinvestigate the Hamiltonian 
of the ionic system, the non-positive-definite problem is 
naturally solved and the phonon Hall conductivity is cal- 
culated exactly by applying the Green-Kubo formula in 
the ballistic thermal transport region. The paper is orga- 
nized as follows. In sec. II we first introduce the model, 
then outline the derivation of second quantization and 
also of thermal conductivity. In sec. HI we present nu- 
merical results some lattices, and also apply our method 
to the sodium chloride (NaCl) crystal in Sec. IV. Finally 
we conclude with a short discussion in sec. V. 



II. THE HAMILTONIAN IN A UNIFORM 
MAGNETIC FIELD 

The Hamiltonian operator of a crystal lattice in the 
presence of a magnetic field is, according to the minimal 
energy principle, given hy^ 



W xx' a/3 

Here Mx and e are the mass and charge respectively of 
the x-th ion in the unit cell, which is assumed to have r 
ions, Ua{lx) is the a-th cartesian coordinate of the dis- 
placement of the a;-th atom in the Z-th unit cell from the 
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equilibrium position Xailx), Pa{lx) is the correspond- 
ing linear momentum operator, (j)cii30',x;l' ,x') are the 
atomic force constants for the crystal and A is the mag- 
netic vector potential. In the presence of Lorentz gauge 
(V • A = 0) the vector potential can be written as a 
linear combintation of equilibrium position of each ion 
seperately Aa{l,x) = '^pCa/3U/3{l,x) where C is an an- 
tisymmetric 3x3 matrix. Pa{l,x) and Ua{l,x) satisfy 
the commutation relation 



[Ua{l,x),Pp{l\x')] ^ ihSapSwSx 



(2) 



Using this commutation relation and substituting the ex- 
pression for A, the Hamiltonian of the system can be 
written in a compact form 



H = ^P^P + ^u^Ku + u^Ap, 



(3) 



where u denotes the column vector of displacements of 
the atoms away from the equilibrium positions for all 
the degrees of freedoms, multiplied by a corresponding 
square root of mass \fM^, p is the associated conjugate 
momenta. Superscript T stands for transpose. The K 
matrix has the following form 

Kax,f}x' [l,x;l ,x) = ^ ^ [A)a^^i3^,dw , (4) 

where the matrix elements of A are the cyclotron fre- 
quency for different masses and is given by 



A 



OLX.j3x' 



(5) 



If we take the magnetic field B along z direction, B ~ Bz, 
then the typical form of A matrix would be 



n 

-n 





(6) 



where =eB/M. The term u^Ap can be treated as on- 
site potential which was introduced phenomenologically 
in all other models to describe PHE but it comes out 
naturally for ionic crystals. Also because of the special 
form of the K matrix {K = cj) — A^) and A being anti- 
symmetric matrix the system is indeed stable for large 
values of magnetic field. 

Using the Heisenberg equation of motion and Bloch's 
theorem we can write the equation of motion for u and 
p in the following matrix form 



Mx = (^X 



where 



A D 

-I A 



(7) 



(8) 



Xcr(k) = (^o-(k),eo-(k))^, Ai<T(k) and eo-(k) are the po- 
larization vectors corresponding to p and u, respectively. 



D(k) = Kii'e'^'^^''~^''> is the dynamical matrix, R; 
is the real-space lattice vector and / is an identity matrix. 

The displacement and momentum operators can be 
written in the following standard second quantization 
form 



ui 



Pi 



k,(T 



2K(k)|iV 



acr(k); 



2|u;.(k)|7V 



a,(k). (9) 



Where k is the wavevector, N is the number of unit cells, 
a is the phonon branch index which runs over both pos- 
itive and negative frequencies, |wo-(k)| — W(j(k)sign(cr) 
and a_o.(k) = aj.(— k). The reason that we have nega- 
tive frequencies is because they are also the eigenvalues 
of Eq.(7). Only both the positive and negative frequency 
eigenmodes form a complete set. Using this complete set 
of eigenvalues and eigenvectors the Hamiltonian can be 
re-expressed as 



7J=i5]ftK(k)|at(k)a,(k) 



(10) 



k,(T 



Using the defination of local energy density and continu- 
ity equation 



(11) 
(12) 



S=^J d'xSix), 
total heat current operator can be written as^"— 

S^ = lY.(^i-^i'>f^i.^'^^'^ (13) 



Ll' 



where ui = pi — Aui. The index a corrensponds to dif- 
ferent cartesian axes. V is the total volume of N unit 
cells. In terms of creation and annihilation operators the 
above expression can be written in an exact form as 



LOa' (k) 



VK(k)|K.(k)| 



,(k)a<^'(k). 



where the F function is defined as 



F;,,(k) = 4(k)^..'(k). 



(14) 



(15) 



is a Hermitian matrix with real diagonal elements 
and is related to the group velocity W(j(k) as F^^{k) — 
2?;o- (k)W(j (k) Off-diagonal elements are in general not 
zero and is responsible for carrying heat in Hall ef- 
fect. The above expression of Sa contains all four possi- 
ble combinations of creation and annihilation operators 
which are ajj(k)aCT'(k), acr(k)a^, (k), ajj(k)a^, (— k) and 
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aCT(— k)acr' (k). The measured heat current J = Tr(pS) 
vanishes in equihbrium. The thermal conductivity tensor 
can be calculated from the Green-kubo formula 



20 



,15.16 



V 
T 



dX 



dt{S''{-iXfi)S''{t))^ 



(16) 



where j3 — l/{kBT), the average is taken over the equi- 
librium ensemble with Hamiltonian H. Substituting the 
expression 5° and using Wick's theorem we obtain 



^ y f;,,(k)xf,v(k)x ("+"^^' 



2,2VT 

niuj') — niuj) 
X— — X 



rj — i{uj — Lo') ' 



(17) 



where n(w) — (e'^'*'^ — 1)^^ is the Bose distribution 
function. We have used the notations lo = a;o-(k), 
u' = tJo-'(k)- r] comes from the damping term e"''* which 
we have added in order to integrate the oscillatory fac- 
tor e^^^~'^ Same phonon branch (cr = a') does not 
contribute to the thermal conductivity. For high tem- 
perature (/3huj < 1) thermal conductivity reaches a con- 
stant value but for low temperature {T — 0) it diverges 
as ±1/T. For a perfect crystal without the magnetic 
field(yl=0) the diagonal elements of thermal conductiv- 
ity Kaa diverge in the form of l/ry , corresponds to infi- 
nite conductivity and the off-diagonal elements Kab are 
zero which is consistent with Fourier's law. For nonzero 
magnetic field (^4 ^ 0) the diagonal elements of k are 
still infinite but off-diagonal elements gives rise to finite 
value of conductivity. The form of Eq.(17) is very similar 
to the anomalous Hall effect for electrons where the sum- 
mation is replaced by an integral of kii. The same form 
was also derived for disordered harmonic crystal and for 
amorphous solidai^ii^. The formula can also be consid- 
ered from a topological point of viewii. The major dif- 
ference of our result with othcrsi^ii^ is in the expression 
for S"^ where we have taken all possible combinations of 
creation and annhilation operators. As a result we need 
to sum over all positive and negative phonon branches. 
We numerically calculate the conductivity and found that 
the effect is almost twice in comparison with Ref. i8|. 



III. CALCULATION ON SOME LATTICES 

Before studying the real ionic-crystal material, we cal- 
culate the PHE on some lattices to get some basic prop- 
erties of the PHE in the ionic lattices. In this section 
we show the calculation of the PHE on the hexagonal 
lattices, triangular lattices, and square lattices. 

In Fig. 1, we give results of dispersion relation for a 
hexagonal lattice with nearest neighbour couplings. The 
primitive vectors for the lattice are (a,0) and (a/2,\/3a/2) 
with a = 1 A. The longitudinal and transverse force con- 
stants are chosen as = 0.144 eV/(uA^) and Kt = 
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FIG. 1: (Color online) (a) Phonon dispersion curve for hexag- 
onal lattice. The angular frequency (both positive and nega- 
tive) of one accoustic mode as a function of kxa with ky — 0. 
The solid curves are for h = and dotted curve are for 
h — 10.0 X 10^^ rad/s. (b) The eigenfrequencies (two accous- 
tic and two optical) as a function of h for fixed wave vector 
ka = (1, 1). 



0.036 eV/(uA^) respectively^. This typical choice is com- 
parable with experimental values. Here h = eB/m which 
has dimension of frequency. In the presence of magnetic 
field the phonon spectrum is shifted (see Fig. 1(a)), and 
the shift is proportional to the applied magnetic field. 
It can also be shown that for lattices having inversion 
symmetry uja{—'k) = ijJaO^) — ~'^-o-(k) even for nonzero 
magnetic field. In Fig. 1(b) we have given a plot both for 
positive and negative eigenfrequencies ujaO^) as a func- 
tion h and we can see that all eigenvalues are real even 
for very large value of h. 

In Fig. 2(a), we show result for the off-diagonal ther- 
mal conductivity K^y as a function of magnetic field for 
three different temperatures T=10K, 50 K, and 100 K. 
It is clear that, at low temperature for small values of 
h, Kxy increases linerly with the magnetic field but for 
large h, it grows slower than linear. For very large h 
it can even goes to negative value. At low temperature 
phonon Hall conductivity can change sign with increas- 
ing magnetic field. This result is consistent with Ref. |^. 
With increasing temperature, the magnitude of negative 
Hall conductivity decreases because when the tempera- 
ture is high, more high energy modes are excited and 
spin-phonon interaction can not easily turn around the 
phonons. 

In Fig. 2(b) we plot the temperature dependence of 
for three different values of h. For small values 
of h phonon Hall effect is almost ncgligablc. But for 
high value of h it increases slowler than linear with 
T and for very high temperature it reaches a constant 
value. For low temperature K~y ^ j^i-ooga. 005,1.012 j-^j, 

h = (10.0, 1.0, 0.1) X 10^2 rad/s respectively(see Fig. 2(c)) 
which is close to 1.0 and consistent with Eq.(17). The di- 
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InT 

FIG. 2: (Color online) Thermal Hall conductivity K^y (a) as a 
function of the magnetic field h (b) as a function of tempera- 
ture T (c) divergence of Ifall conductivity at low temperature 
for a hexagonal lattice. 



vergence of n^y as 1 /T for T — >■ is not unphysical as it 
shows the ballistic behaviour for which the longitudinal 
conductivity Kxx is always infinite at any temperature. 
However if we measure current it will be zero at T — 0. 
And from another point of view, when T — > 0, the value 
of KxyT is constant; because of AT < T in linear response 
region, the measurable quantity of transverse heat cur- 
rent — —KyxAT/L ( we know Kyx — —Hxy ) will 
not be divergent and tend to zero for the periodic sys- 
tem L — > oo. Furthermore we know that for the ballistic 
thermal transport, because of the zero temperature gra- 
dient in the system, the thermal conductivity is infinite 
while the thermal conductance is finite and meaningful. 
Therefore, the 1/T divergence of Kxy indeed embodies 
the ballistic picture of PHE at very low temperature. 

In Fig 3. we show results for a triangular lattice. We 
have taken the same values of the parameters for the cal- 
culation. In this case also the off-diagonal thermal con- 
ductivity shows similar behaviour as in hexagonal lattice. 
At very low temperature thermal conductivity changes 
sign. The temperature depedence is also similar. At low 
temperature conductivity goes as —1/T. 

For a square lattice with only the nearest neighbour 
coupling there is no PHE. Generally the off-diagonal term 




h(10'^rad/s) T (K) 

FIG. 3: (Color online) Thermal Hall conductivity K^y (a) as 
a function of the magnetic field h and (b) as a function of 
temperature T for a triangular lattice. 

of K will be zero if the reflection symmetry is not bro- 
ken. If there exists a symmetry operation S such that 
SKS-^ = K and SAS-^ = -A then Hah = for a^h. 
In this case the spring constant matrix K is diagonal and 
is invariant under reflection in a; or y direction, hence it 
satisfies the above relation. If we consider next-neighbour 
coupling then the matrix K will not have mirror reflec- 
tion symmetry and the PHE can be observed which is 
consistent with Ref. 



IV. THE PHE ON THE SODIUM CHLORIDE 
CRYSTAL 

After the properties of the PHE are well studied for 
the model lattices, we turn to the calculation for the real 
ionic crystal NaCl. NaCl crystal can be described as fee 
lattice with primitive vectors ai — |-(y-fz), ~ f (x+z) 
and as = ^ (x + y) and two point basis consisting of a 
sodium ion at and a chlorine ion at the center of the 
conventional cubic cell, |(x-f y-l-z). Here a is the lattice 
constant which is 5.63 A for NaCl. 

We calculate the force constant matrix using "General 
Utility Lattice Program" (GULP)^ by considering a unit 
cell. The interatomic potential we used includes both 
Coulombic and non-Coulombic terms 

yij{fi,rj) — ^— — h Aijexp{—rij / Bij) g^, 

(18) 

where rij is the interatomic distance between ions i and 
j and Zi is the effective charge of the ith ion. The non- 
Coulombic part of the potential is known as Buckingham 
potential. Aij and Bij are parameters for the repulsive 
interaction and Cij is van der Waals constant. The val- 
ues of the parameters are given in Table 1?^ . Cutoff 
distance 12 A is used in the calculation. With this po- 
tential model we went upto the third nearest neighbours 
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TABLE I: Parameters of the potential of NaCl 



Source 


A (ev) 


B(A) 


C (e«)A6 


Na-Na 


487 


0.23768 


1.05 


Na-Cl 


14514 


0.23768 


6.99 


Cl-Cl 


405774 


0.23768 


72.40 




FIG. 4: (Color online)Thermal Hall conductivity Kxy (a) as 
a function of the magnetic field B (b) and as a function of 
temperature T for NaCl. 



and cut the force constant matrix properly in order to 
get the dynamical matrix. The choice of this type of unit 
cell gives us dynamical matrix of size 6x6. 

In Fig. 4(a), we give the off-diagonal thermal conduc- 
tivity Kxy as a function of B for temperature T = 100 -ftT 
and T — 50 K. We note that the sign of Kxy is nega- 
tive. The growth of Kxy in this case, for small values of 
magnetic field is very fast, then it reaches a maximum 
and then decreases to zero with increasing magnetic field 
and even goes to positive value for high value of magnetic 
field. 

In Fig. 4(b), we display the temperature dependence 
of the off-diagonal thermal conductivity Kxy for B = 0.3 
and B = 0.5 MTesla. Upto r=40K the off diagonal 
thermal conductivity is close to zero, meaning there is 
no heat flow. When T — it also diverge like l/T. But 
it gradually increases and finally saturates at high tem- 
perature. In all our calculations we have used rj of order 
10~^Kl to get rid of the sigularity due to the accoustic 
phonons at F point. 

From the above calculation in the nonmagnetic ionic 
crystal NaCl, we observe the PHE although there is no 
magnetization in the system. The existence of PHE is 
due to the nonzero term of A in Eq. ^ . Both the positive 
and negative ions in the crystal contribute to A, and they 
can not calcel each other. Since in the term A, the charge 
to mass ratio (e/m) is very small for nonmagnetic ionic 
crystals NaCl, a very strong magnetic field is required 



to observe PHE in this case. However in the paramag- 
netic crystal the contribution from electrons dominates 
the term A, which increases about 10^ times and thus the 
PHE can be comparable to the experimental data. For 
ferromagnetic insulating ionic crystal, it will be increased 
much more because of the large magnetization which will 
contribute to large spin-phonon interaction. 

V. CONCLUSION 

A theory for phonon Hall effect in general ionic crystals 
based on a ballistic lattice dynamic model is developed. 
An exact formula of phonon Hall conductivity formula 
from Green-Kubo approach is derived. It is shown that 
the effect can be present in the presence of a very strong 
magnetic field for general ionic crystals without magne- 
tization. Due to the tiny charge-to-mass ratio, the PHE 
is so small that it is difficult to be measured experimen- 
tally. However the theoretically confirmed existence of 
PHE can give us deep understanding of the PHE and 
guide us the direction to do the experiments. The ferro- 
magnetic insulator is a good candidate for the PHE. We 
hope some artificial lattices such as cool atoms in opti- 
cal trap can generate similar kind of effect. We believe 
that the result for Kxy is robust against nonlinear pertur- 
bation, but it is interesting to work out a theory where 
nonlinearity is taken into account. 
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Appendix A: Second Quantization 

For the negative phonon branches we define the polar- 
ization vectors as following 

6_^(-k) = 6:(k), (Ai) 

L^_,(~k) = -c^,(k). (A2) 

This defination is consistent with Eq. (7). Also the cre- 
ation and annihilation operators are defined as for ct > 
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a_,(k) = aii-k), (A3) 
and the time dependence is given by for cr > : 



(A4) 
(A5) 



The commutaion relation between the operators is 

[a„{k),a,,{k')]=5{k + k')S{a + a')sign{a). (A6) 

The matrix A4 in Eq.(8) is not a normal matrix. 
Hence we need to consider both the left and right 
eigenvectors. However it can be easily shown that 
they are not really independent. For a given value 
of uja{k) and corresponding right eigenvector Xcr(k) = 
(/LtCT(k), ecr(k))-^ one can choose the left eigenvector 
as^X^k) = (et(k)/(-2*^.(k)),-A^Uk)/(-2*L.,(k))). 
With this choice of eigenvectors the normalization can 
be done according to 



4(k)-e.(k) 



-et(k)-i-e.(k) = l. (A7) 



Wo-(k) 

The Hamiltonian in Eq. (3) can be written as 



LI' 



' I -IS, I, A6u' / ^ ' 



where 

Xi = 

X, = I 



2K(k)|7V 
^(-2zc.,(k))x.(k)e-^^'-'' > 



a<T(k); 



2K(k)|iV 



4(k), 



(A9) 



where we have used the definition for ui and pi and 
also the Hermitian property of the operator. Now us- 
ing the definition of dynamical matrix, the identity that 
Y,i e*('''-'^)-^' = NS{k' - k), and the orthogonality rela- 
tion between left and right eigenvector x^{k) ■ Xo-'(k) = 
Sa-a', the Hamiltonian reduces to the form of Eq. (11). 
The canonical commutation relation can be shown to be 
satisfied using the completeness of x, as 



-ihSwI 



(AlO) 



Appendix B: Positive definiteness of the 
Hamiltonian 

To prove that the Hamiltonian is positive definite we 
can write the Hamiltonian as 



H 





(Bl) 



K 2A 

with K = (j) — A^. We assume that (j) is positive definite. 
Since transpose of a positive definite matrix is positive 
definite we can as well consider the matrix 



(B2) 



which is a symmetric matrix. Hence it can be diagonal- 
ized by similarity transformation and can be written as 

Q = SVS'^ where V =diag(Ai, Aa, ). Now if A; > 

V i then 2? 2 exists and we can write Q — B^B with 
B = D'^ S'^ . So a real symmetric matrix has non- negative 
eigenvalues ifi^ it can be factorized as B^B. The above 
matrix Q can be factorized according to 




B 



-A I 



^K + A^Y ) 
Hence the Hamiltonian is positive definite. 
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